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Abstract 

We report measurements of the first and second moments of the hadronic invariant mass squared 
distribution, (M^) and {{M\ - (M^)) 2 ), in B -> X c iu decays for minimum lepton momenta 
ranging from 0.7 to 1.5 GeV/c in the B meson rest frame. The measurement uses BB events in 
which the hadronic decay of one B meson is fully reconstructed and the semileptonic decay of the 
other B is inferred from the presence of an identified lepton. These results are obtained from a 
140 fb _1 data sample collected near the T(45) resonance with the Belle detector at the KEKB 
asymmetric energy e + e~ collider. 
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Recently, there have been intense theoretical and experimental efforts toward predict- 
ing [l[ HI HI and measuring [3, 13, 0] the moments of the hadronic invariant mass squared 
distribution in B — > X c iv decays. The idea is that using the Operator Product Expansion 
(OPE) [3], the hadronic mass moments (and other inclusive observables in B decays) can 
be predicted in terms of the 6-quark mass m;, and other non-perturbative parameters. Con- 
versely, by measuring the moments of B decay spectra and the semileptonic B decay rate, 
one can then extract these parameters and \Vd\ j9[. 

The analysis is based on the data recorded with the Belle detector [13] at the asymmetric 



energy e + e~ collider KEKB jl ljj . operating at a center-of-mass (cm.) energy near the 
T(45) resonance. KEKB consists of a low energy ring (LER) of 3.5 GeV positrons and 
a high energy ring (HER) of 8 GeV electrons. The Belle detector is a large-solid-angle 
magnetic spectrometer consisting of a three-layer silicon vertex detector (SVD), a 50-layer 
central drift chamber (CDC), an array of aerogel threshold Cerenkov counters (ACC), a 
barrel- like arrangement of time-of-flight scintillation counters (TOF), and an electromagnetic 
calorimeter comprised of CsI(Tl) crystals (ECL) located inside a super- conducting solenoid 
coil that provides a 1.5 T magnetic field. The responses of the ECL, CDC (dE/dx) and 
ACC detectors are combined to provide clean electron identification. Muons are identified 
in the instrumented iron flux-return (KLM) located outside of the coil. Charged hadron 
identification relies on the information from the CDC, ACC and TOF sub-detectors. 

The T(4S) dataset used for this study corresponds to an integrated luminosity of 140 fb _1 , 
or about 152 million BB events. Another 15 fb _1 taken 60 MeV below the resonance are 
used to subtract the non-BB (continuum) background. Full detector simulation based on 
GEANT [13] is applied to Monte Carlo (MC) simulated events. The size of the MC samples 
is equivalent to about 2.4 times the integrated luminosity. At the generator level, the decay 
B -> D*iu is simulated using a HQET-based model Q. The ISGW2 model Q is used 
for the decays B — ► Div and B — > D**£u. The non-resonant B — > D^ixtv component 
is generated according to the model of Goity and Roberts [15]. QED bremsstrahlung in 
semileptonic decays is simulated by the PHOTOS package [16j. 



After selecting hadronic events [17|, we fully reconstruct the hadronic decay of one 
B meson (5 tag ) using the decay modes B + -> Z)(* )0 7r + , £>W°p+, D^°af and B° -> 
DW-7T+ D^~p+, D<d-<4 [18]. Pairs of photons satisfying £ 7 > 50 MeV and 117 MeV/c 2 < 
771(77) < 150 MeV/c 2 are combined to form 7r° candidates. K s mesons are reconstructed 
from pairs of oppositely charged tracks with invariant mass within ±30 MeV/c 2 of the 
Kg mass and decay vertex displaced from the interaction point. Candidate p + and p° mesons 
are reconstructed in the ir + n° and 7t + tt~ decay modes, requiring their invariant masses to 
be within ±150 MeV/c 2 of the p mass. Candidate mesons are obtained by combin- 
ing a p° candidate with a charged pion and requiring an invariant mass between 1.0 and 
1.6 GeV/c 2 . D° candidates are searched for in the K~ir + , K~7t + tt°, K~tt + 7t + tt~ , K s tt + 7t~ 
and KgTi decay modes. The K^n + n + and K s n + modes are used to reconstruct D + mesons. 
Charmed mesons are selected in a window corresponding to ±3 times the mass resolution 
in the respective decay mode. D* + mesons are reconstructed by pairing a charmed meson 
with a low momentum pion, D* + — > D°tt + , D + tt°. The decay modes D*° — > D°tt and 
D*° — > D 0/ ~f are used to search for neutral charmed vector mesons. 

For each B tag candidate, the beam-constrained mass M^ c and the energy difference AE 
are calculated 



M bc = V^beam) 2 - {p B ) 2 , AE = E B - E heam , (1) 

where i?beam is the beam energy in the cm. system and p B and Eg are the 3-momentum and 
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the energy of the B tag candidate in the same frame, respectively. The signal region for B tSLg 
is defined by the selections Mb c > 5.27 GeV/c 2 and \AE\ < 50 MeV. If multiple candidates 
are found in a single event, the best candidate is chosen based on AE and other variables. 

Semileptonic decays of the other B meson (-B S i g nai) are selected by searching for an iden- 
tified charged lepton (electron or muon) within the remaining particles in the event. Events 
with multiple identified leptons are rejected. In the lepton momentum range relevant to this 
analysis, electrons (muons) are selected with an efficiency of 92% (89%) and the probability 
to misidentify a pion as an electron (a muon) is 0.25% (1.4%) [19|,|20j]. For B + tags, we require 
that the lepton has a charge sign compatible with a prompt semileptonic decay (Qi-Qb < 0, 
where Qe and Qb are the charges of the lepton and of -Stag, respectively). In electron events, 
we attempt to recover bremsstrahlung photons by searching for a E 1 < 1 GeV photon within 
a 5° cone with respect to the electron direction. If such a photon is found, it is merged with 
the reconstructed electron candidate. 

Charged and neutral particles in the event associated neither with B tSug nor with the 
charged lepton are assigned to the hadronic X system. The missing 4-momentum in the 
event is calculated, assigning the pion mass to all charged particles except identified kaons, 

Pmiss = G»LER + PHER) - PB tag ~Pl~PX , (2) 

where the indices LER and HER refer to the colliding beams. As only the neutrino in 
B — > X c iv should be missing in the event, the missing mass is required to be consistent 
with zero, |M^ iss | < 3 GeV 2 /c 4 . To improve the resolution in M x , we constrain the neutrino 
mass to zero, p u = (\p m i S s\,Pmiss), an d recalculate the 4-momentum of the X system, 

Px = (Pler + Pher) - PB t , g -Pi-Pv ■ (3) 

The M\ resolution (defined as half width at the half maximum) obtained in this way is 
about 800 MeV 2 /c 4 . 

For the rest of the analysis, the remaining events are divided into four sub-samples, de- 
pending on the charge of B tag (B + , B°) and on the lepton type (electron, muon). In each 
of these sub-samples and for each lepton momentum threshold considered in the analysis 
(pi > 0.7, 0.9, 1.1, 1.3 and 1.5 GeV/c 0), the backgrounds in the M\ distribution are 
determined, taking into account contributions from the following sources: continuum back- 
ground, BB events with a misreconstructed B tag candidate, and background from secondary 
or fake leptons. The background shapes in M\ are determined from the MC simulation, 
except for the continuum where the off-resonance data is used. The shape of the fake muon 
background is corrected by the ratio of the pion fake rate in the experimental data over the 
same quantity in the MC simulation, as measured using K% — > 7r + 7r~ decays. The contin- 
uum background is scaled by the on- to off-resonance luminosity ratio, taking into account 
the cross-section difference. The combinatorial B tSig background is normalized using the 
5.20 GeV/c 2 < M^ c < 5.25 GeV/c 2 sideband. The normalization of the secondary or fake 
lepton background is found from the real data by fitting the lepton momentum distribution. 
The purity of the B — > X c iv signal depends on the sub-sample and the lepton momentum 
threshold, typical values being around 75%. Table shows the numbers of signal events and 
purities for each combination of B tag charge, lepton type and lepton momentum threshold. 

In each of the four sub-samples and for each lepton momentum threshold, the M\ distri- 
bution is measured in 39 bins in the range from to 13 GeV 2 /c 4 (bin width 0.333 GeV 2 /c 4 ) 
and, after subtraction of all backgrounds, the finite detector resolution in M\ is unfolded 
using the Singular Value Decomposition (SVD) algorithm ji^], as illustrated in Fig. [T] The 
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p^m B + electron B + muon B electron B muon 

0.7 3893 ± 79 (72.3%) 3626 ± 84 (68.7%) 2212 ± 64 (66.4%) 2154 ± 64 (65.u/o ) 

0.9 3659 ± 75 (74.0%) 3484 ± 78 (71.1%) 2072 ± 58 (66.4%) 2067 ± 58 (70.2%) 

1.1 3285 ± 70 (75.2%) 3159 ± 73 (72.9%) 1886 ± 53 (66.4%) 1925 ± 52 (76.1%) 

1.3 2742 ± 64 (75.9%) 2740 ± 66 (75.0%) 1595 ± 46 (66.4%) 1632 ± 47 (79.4%) 

1.5 2152 ± 56 (77.7%) 2132 ± 56 (76.8%) 1195 ± 39 (66.4%) 1297 ± 41 (83.0%) 



TABLE I: Number of B — > X c lv signal and signal purity in the four sub-samples, as a function 
of the lepton momentum threshold. The yields are quoted with their statistical uncertainty; the 
corresponding signal purity is given in brackets. 
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FIG. 1: (a) Measured and (b) unfolded M\ distribution for <p\ > 0.7 GeV/c. On the left plot, 
the continuum-subtracted M\ distribution is shown by points with error bars. The hatched his- 
togram corresponds to the background from secondary or fake leptons; BB events in which Stag 
is misreconstructed are shown by the double-hatched histogram. The right plot is the result of the 
unfolding. In both plots, contributions from B + and B° tags, and from electron and muon events 
are added. 

unfolded distribution has 15 bins in the range from M|, to about 15 GeV 2 /c 4 (bin width 
1 GeV 2 /c 4 , except around the narrow states D, D*, D\ and D^). We calculate the first and 
second moment, (M|) and ((M| — (M x )) 2 ), for each unfolded M\ spectrum separately, 
after applying a small correction for different bin-to-bin efficiencies. The final results for 
a given lepton momentum threshold are obtained by taking the average over the four sub- 
samples. The unfolding, moment calculation and averaging procedure has been studied on 
MC simulated events and no significant bias has been found. 

The results for the first and second hadronic mass moment are shown in Fig. |21 and 
Tables [H] and IIHI All results are preliminary. Note that the moment measurements for 
different lepton momentum thresholds are highly correlated due to overlapping data samples. 
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FIG. 2: (a) First and (b) second hadronic mass moment, (M^) and ((M| - (M|)) 2 ), for different 
lepton threshold momenta. The error bars indicate the statistical and total errors. Note that the 
individual moments are highly correlated. All results are preliminary. 





{Ml) 


detector/ 


unfolding 


X c moc 


(GeV/c) 


(GeV 2 /c 4 ) 


background 






0.7 


4.383 ±0.037 ±0.051 


0.047 


0.004 


0.020 


0.9 


4.330 ± 0.033 ± 0.041 


0.036 


0.009 


0.018 


1.1 


4.277 ±0.029 ±0.035 


0.028 


0.014 


0.016 


1.3 


4.173 ±0.028 ±0.037 


0.026 


0.021 


0.016 


1.5 


4.132 ±0.030 ±0.031 


0.019 


0.018 


0.016 



TABLE II: First hadronic mass moment (Mj^) for different lepton threshold momenta. The first 
error on (M^) is statistical and the second is the estimated systematic uncertainty. The right-most 
three columns show the different components of the systematic error. All results are preliminary. 

We have estimated the correlations due to this overlap in Table ITVl 

We consider three sources of systematic error, shown separately in columns three to five 
of Tables HU and IIHI the uncertainty related to the detector modeling and the background 
subtraction, the uncertainty related to unfolding and the moment extraction procedure, and 
the uncertainty related to the X c model in the MC simulation. The total systematic error 
is the quadratic sum of these three components. 

The uncertainty related to the detector modeling and the background subtraction is 
estimated by varying the normalizations of the different background components and some 
of the selections used in the analysis (namely the B taig signal region and the requirement on 
M^ iss ). The uncertainty related to unfolding and moment extraction is obtained by varying 
the effective rank parameter in the SVD algorithm, dis- and enabling the bin-to-bin efficiency 
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p* mhl {{Ml - (M 2 )) 2 ) 
(GeV/c) (GeV 4 /c 8 ) 

0.7 1.605 ±0.179 ±0.322 

0.9 1.317 ±0.144 ±0.236 

1.1 1.013 ±0.113 ±0.161 

1.3 0.634 ±0.095 ±0.105 

1.5 0.591 ±0.110 ±0.088 



detector/ unfolding X c model 
background 

0.233 0.184 0.123 

0.174 0.116 0.110 

0.125 0.047 0.089 

0.079 0.028 0.064 

0.035 0.068 0.043 



TABLE III: Same as Table ITU for the second hadronic mass moment {{Ml 



(Mi)) 2 }. 



Pmin 

(GeV/c) 



{Ml) 



0.7 0.9 1.1 1.3 1.5 



mi 

0.7 0.9 



- (Mi)) 2 ) 
1.1 1.3 



1.5 



{Ml) 



0.7 
0.9 
1.1 
1.3 
1.5 



1.000 0.922 0.807 0.620 0.533 
1.000 0.875 0.672 0.579 
1.000 0.768 0.661 
1.000 0.861 
1.000 



0.782 0.720 
0.705 0.781 
0.561 0.622 
0.347 0.384 
0.342 0.379 



0.622 
0.674 
0.770 
0.476 
0.470 



0.471 
0.511 
0.583 
0.759 
0.749 



0.415 
0.450 
0.514 
0.670 
0.778 



0.7 
0.9 

{{Ml - {Ml)f) 1.1 
1.3 
1.5 



1.000 0.902 
1.000 



0.728 
0.807 
1.000 



0.456 
0.506 
0.627 
1.000 



0.440 
0.487 
0.603 
0.963 
1.000 



TABLE IV: Estimated correlation coefficients for the different moment measurements due to over- 
lapping data samples. 

correction and changing the binning of the unfolded distribution. The X c model uncertainty 
is determined by varying the fractions of B -> D*£u, B -»■ Div and B -> D**/D^tt£u within 
±10%, ±10% and ±30%, respectively, and summing the individual variations in quadrature. 
These ranges of variation roughly correspond to the experimental uncertainties [23| in the 
isospin averaged branching ratios. 

In summary, we have measured the first and second moments of the hadronic invariant 
mass squared distribution, {Ml) and {{Ml — {Ml)) 2 ), in B — > X c lv decays for minimum 
lepton momenta ranging from 0.7 to 1.5 GeV/c in the B rest frame. The results obtained 
are compatible with theoretical expectations |l| and recent measurements from other exper- 
iments p, 0| • In addition, we have estimated the correlations of the moment measurements. 
These measurements can be used as input to a global fit analysis, which is expected to lead 
to an improved determination of the CKM matrix element \V c b\. 
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